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In the last time the cold dark matter (CDM) model has suggested more and more that it is not 
able to describe all the properties of nearby galaxies that can be observed in great detail as well 
as that it has some problems in the mechanism by which matter is more rapidly gathered into 
large-scale structure such as galaxies and clusters of galaxies. In this work we revisit an alternative 
model, the scalar field dark matter (SFDM) model, which proposes that the galactic haloes form 
by condensation of a scalar field (SF) very early in the Universe, i.e., in this model the haloes of 
galaxies are astronomical Bose-Einstein Condensate drops of SF. On the other hand, large-scale 
structures like clusters or superclusters of galaxies form similar to the ACDM model, by hierarchy, 
thus all the predictions of the ACDM model at cosmological scales are reproduced by SFDM. This 
model predicts that all galaxy haloes must be very similar and exist for higher redshifts than in the 
ACDM model. In the first part of this review we revisit the cosmological evolution of SFDM model 
with a scalar potential m 2< I> 2 /2 + A$ 4 /4 with two different frameworks: the field and fluid approach. 
We derive the evolution equations of the SF in the linear regime of perturbations as well. The scalar 
fluctuations have an oscillating growing mode and therefore, this kind of dark matter could lead 
to the early formation of gravitational structures in the Universe. We also revisit how BEC dark 
matter haloes exhibit a natural cut of the mass power spectrum. In the last part, we study the 
core central density profiles of BEC dark matter haloes and fit high-resolution rotation curves, we 
show a sample of some low surface brightness galaxies. The mean value of the logarithmic inner 
density slopes is a = —0.27 ± 0.18. Using a model independent new definition of the core in the 
BEC density profile, we show that the recent observation of the constant dark matter central surface 
density can be reproduced. We conclude that in light of the difficulties that the ACDM model is 
currently facing the SFDM model can be a worthy alternative to keep exploring further. 

PACS numbers: 95.36. +x,03.65.-w 



INTRODUCTION 

Nowadays the most accepted model in cosmology 
which explains the evolution of the Universe is known 
as ACDM. In this model 4 per cent of the total content 
of the Universe is baryonic matter, 22 per cent is non- 
baryonic dark matter (DM) and the rest is in some form 
of cosmological constant. ACDM has achieved several 
observations with outstanding success, like for example 
the fact that the cosmic microwave background radia- 
tion can be explained in great detail and that it provides 
a framework within one can understand the large-scale 
isotropy of the Universe and important characteristics on 
the origin, nature and evolution of density fluctuations 
which are believed to give rise to galaxies and other cos- 
mic structures. There remain, however, certain conflicts 
at galactic scales, like the cusp profile of central densi- 
ties in galactic halos, the overpopulation of substructures 
predicted by A-body numerical simulations which are an 
order of magnitude larger than what has been observed, 
etc (see for example [TSJ, [50], [75], and [ST]). 

There are other important issues of the ACDM model 
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which require further considerations. Observations point 
out to a better understanding of the theory beginning 
with the Local Void, which contains far fewer galaxies 
than the expected. Another problem arises for the so 
called pure disk galaxies, which do not appear in numer- 
ical simulations of structure formation in the standard 
model. These problems would be solved, if the structure 
grew faster than it does in the standard paradigm [80] . 

On the other hand, Lee & Komatzu [59] also found that 
the collision velocity of 3000km/s at R200 for the Bullet 
Cluster is very unlikely and might no be explained by 
the standard model. A final example of inconsistencies 
can be seen in a paper by Thomas, Abdalla & Lahav 
[95], who found anomalies in the mass power spectrum 
obtained by SDSS and the one obtained with the CDM 
model. With these and other results it seems necessary 
to change the ACDM paradigm to try and explain the 
formation of structure in the Universe. 

The incorporation of a new kind of DM, different from 
the one proposed by the ACDM model into the big bang 
theory holds out the possibility of resolving some of these 
issues. Several authors have proposed interesting alter- 
natives in where they try to solve the latter explained 
difficulties that the CDM scenario seem not solve. Some 
of these models consider warm dark matter ([19] ■ |101j ). 
self interacting dark matter ([61], [93]), and other exotic 
scenarios as braneworld models 32 . In fact, there are 
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models that do not include DM but instead modify the 
Newtonian force law [74, 90J . Other alternative scenario 
that has received much attention in the last years is 
the scalar field dark matter (SFDM) model. The main 
idea is simple (see Guzman & Matos [35]), the nature 
of the DM is completely determined by a fundamental 
scalar field The SFDM model proposes that galac- 
tic haloes form by Bose-Einstein condensation of a scalar 
field (SF) whose boson has an ultra-light mass of the or- 
der of m ~ lCU 22 eV. From this mass it follows that the 
critical temperature of condensation T c ~ 1/to 5 / 3 ~TeV 
is very high, therefore, they form Bose-Einstein Conden- 
sates (BEC) drops very early in the Universe. In ad- 
dition, the Compton length A c = 2nh/m associated to 
this boson is about ~ kpc that corresponds to the dark 
halo-size of typical galaxies in the Universe. Thus, it has 
been proposed that these drops are the haloes of galax- 
ies (see Matos & Urena [70] ). i.e., that haloes are huge 
drops of SF. On the other hand, big structures form like 
in the ACDM model, by hierarchy [70] [03], thus, all suc- 
cessful predictions of the standard model at large scales 
are well reproduced by SFDM. In other words, in the 
SFDM model the haloes of galaxies do not form hierar- 
chically, they are formed at the same time and in the 
same way when the Universe reaches the critical temper- 
ature of condensation of the SF. From this it follows that 
all galaxies must be very similar because they formed in 
the same manner and at the same moment [58) . There- 
fore, from this paradigm we have to expect that there 
exist well formed galaxy haloes at higher redshifts than 
in the ACDM model. In this model the scalar particles 
with that ultra light mass are such that their wave prop- 
erties avoid the cusp problem and reduce the high number 
of small satellites by the quantum uncertainty principle 
[46]. [63] which is another problem that is still present in 
the ACDM model |33J [5U1 [S5J? ]■ Summarizing, it is 
remarkable that with only one free parameter, the ultra- 
light scalar field mass (to ~ lCT 22 eV), the SFDM model 
fits: 

1. The evolution of the cosmological densities [7Tj . 

2. The acoustic peaks of the cosmic microwave back- 
ground [55] . 

3. The scalar field has a natural cut off, thus the sub- 
structures in clusters of galaxies is avoided natu- 
rally. With a scalar field mass of ~ 10 _22 eV 
the amount of substructure is compatible with the 
observed one [41) ] 170 ] IM] . 

4. We expect that SFDM forms galaxies earlier than 
the CDM model, because they form BECs at a crit- 
ical temperature T c >>TeV. So if SFDM is right, 
we have to see big galaxies at high redshifts with 
similar features [7T] [03] . 



5. The rotation curves of big galaxies and LSB galax- 
ies [4] H3J [42] [60] [85] . 

6. With this mass, the critical mass of collapse for a 
real scalar field is just 1O 12 M , i.e., the one ob- 
served in galaxies haloes pQ. 

7. The observed properties of dwarf galaxies, i.e., the 
minimum length scale, the minimum mass scale, 
and their independence from the brightness [58] . 

8. And recently it has been demonstrated that if 
the scalar mass is m ~ 10~ 22 eV, the SFDM 
haloes would have cores large enough to explain 
the longevity of the cold clump in Ursa Minor and 
the wide distribution of globular clusters in Fornax 
(see Lora et al. [53]). 

The idea was first considered by Ji S. U. & Sin S. J. 
[37], Sin [01] and independently by Guzman & Matos 
[371 135] , Lee & Koh [55J , Matos & Guzman [HJJ suggested 
bosonic dark matter as a model for galactic halos [28] . 
In the BEC model, DM haloes can be described, in the 
non-relativistic regime, as Newtonian gravitational con- 
densates made up of ultra light bosons condensated in 
a single macroscopic wave function. Thus, BECs haloes 
can be described as a coherent scalar field $. 

Several authors have introduced a dynamic scalar field 
with a certain potential U(<E>) as a candidate to dark 
matter, although there is not yet an agreement for the 
correct form of the potential of the field. One interesting 
work pointing this way was done by [68] and indepen- 
dently by 88. where they used a potential of the form 
V($) = V [cosh — 1] to explain the core density 
problem for disc galaxy halos in the ACDM model (see 
also PH], [53] and [31]) ■ Guzman & Matos [3SJ presented 
a model for the DM in spiral galaxies, in which they sup- 
posed that DM is an arbitrary SF endowed with a scalar 
potential. Other scalar potential widely used to describe 
DM is V($) = to 2 $ 2 /2 [HJ El]. This potential is very 
interesting because it can mimic the cosmological evolu- 
tion of the Universe predicted by the ACDM model. If 
we consider a SF self-interaction, we need to add a quar- 
tic term to the SF potential EHZl HI] , in this case the 
equation of state of the SF is that of a polytope of index 
n=l (see [3J [33J M M 103] ) • 

Different issues of the cosmological behavior of the 
SFDM/BEC model have been studied (see for example 

[2 E EH 021 Ell US] S01 ESS Mi EOHza [86] ia E31 ESJ 

11001 1103) 1 . Hu, Barkana & Gruzinov [30] proposed fuzzy 
dark matter that is composed of ultra light scalar parti- 
cles who are initially in the form of a BEC. In their work 
Woo & Chiueh [103] used a bosonic dark matter model to 
explain the structure formation via high-resolution simu- 
lations, Ureha-Lopez 98, 99] reviewed the key properties 
that may arise from the bosonic nature of SFDM models. 
Recently Harko [33] [33] developed a further analysis of 
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the cosmological dynamics of SFDM/BEC as well as of 
the evolution of their fluctuations (see also [24]). In the 
same direction, Suarez & Matos [94] studied, the growth 
of scalar fluctuations and large-scale structure formation 
with the fluid approach of the SFDM/BEC model. 

On the other hand, SFDM/BEC model has provided 
to be a good candidate for DM haloes of galaxies in the 
Universe because it can explain many aspects where the 
standard model fails ( [H OH S3 GSMl H IBS] ) . In ad- 
dition, many numerical simulations have been performed 
to study the gravitational collapse of SFDM/BEC model 
[13 EE3 HOI 1351 13SHH] ■ Chavanis [H [53] found approxi- 
mate analytical expression and numerical solutions of the 
mass-radius relation of SFDM/BEC haloes. Recently, 
Rindler-Daller & Shapiro [53] give constraints on the bo- 
son mass to form and mantain more than one vortex 
in SFDM/BEC haloes. These constraints are in agree- 
ment with the ultra light mass founded in previous works 
(see also [1H1 1104] ). Lately, Lora et al. [62 performed 
N-body simulations to study the dynamics observed in 
Ursa Minor dwarf galaxy. They model the dark matter 
halo of Ursa Minor as a SFDM/BEC halo to establish 
constraints for the boson mass. Moreover, they intro- 
duce a dynamical friction analysis with the SFDM/BEC 
model to study the wide distribution of globular clusters 
in Fornax. An overall good agreement is found for the 
ultra light mass ~ 10 _22 eV of bosonic dark matter. 

The main objective of this review is to introduce the 
framework of the SFDM/BEC model that assumes that 
DM is a scalar field that involves a self-interacting po- 
tential of the form U($) = m 2 $ 2 /2 + A$ 4 /4, where 
m ~ 10 _22 eV is the mass of the scalar field, [55], [70] 
and 06]. 

We solve the Friedmann equations for the SFDM/BEC 
model to show that it behaves just like the CDM model. 
Also, we show that the SFDM/BEC predicts galaxy for- 
mation earlier than the CDM model, because they form 
BEC at a critical temperature T c » TeV. So, if SFDM 
is right, this would imply that we have to see big galaxies 
at high redshifts. In order to do this, we study the den- 
sity fluctuations of the scalar field from a hydrodynamical 
point of view, this will give us some information about the 
energy density of dark matter halos necessary to obtain 
the observational results of large-scale structure. Here 
we will give some tools that might be necessary for the 
study of the early formation of structure. 

In section [] we study the cosmological evolution of the 
SF with a field and a fluid approach. We obtain a numer- 
ical solution for the scalar density as well as analytical 
expressions for the kinetic and potential energies of the 
SF. Also, we study the evolution for the density contrast 
of scalar fluctuations in the linear regime. In section |we 
study the rotation curves and the cusp/core discrepancy 
with SFDM haloes to compare with the predictions of 
CDM. Finally we give our conclusions in section]] 



THE SCALAR FIELD DARK MATTER 
COSMOLOGY 

The framework of the standard cosmological model is 
of a homogeneous and isotropic Universe whose evolution 
is best described by Friedmann's equations that come 
from general relativity and whose main ingredients can be 
described by fluids whose characteristics are very similar 
to those we see in our Universe. Of course, the Universe is 
not exactly homogeneous and isotropic but this standard 
model does give us a framework within which we can 
study the evolution of the expansion rate of the Universe 
as well as the evolution of small fluctuations that give rise 
to gravitational structures in the Universe as galaxies and 
clusters of galaxies. 

In this section, we study the cosmological dynamics of 
the SFDM/BEC model and the evolution of their scalar 
perturbations with two different frameworks: the field 
and fluid approach. 



Background Universe: The Field Approach 

We use the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric with scale factor a(t). Our background 
Universe is composed by SFDM ($o(£)) endowed with a 
scalar potential V = U(<I>o), radiation (z), neutrinos (v), 
baryons (&), and a cosmological constant (A) as dark en- 
ergy. We begin by recalling the basic background equa- 
tions. From the energy-momentum tensor T for a scalar 
field, the scalar energy density Tg and the scalar pressure 
Tj are given by 

rf = -p*o=-(l&a + v), (1) 

Ti = P* =(~& -V\s i j , (2) 

where the dots stand for the derivative with respect to the 
cosmological time and Sj is the Kronecker delta. Thus, 
the Equation of State (EoS) for the scalar field is P<j> = 
w$ p$ with 

m + v 

The radiation fields, the baryonic component and the 
cosmological constant are represented by perfect fluids 
with baryotropic equation of state P 7 = (7 — l)p 7 , where 
7 is a constant, < 7 < 2. For example, ^ z — j v = 4/3 
for radiation and neutrinos, 75 = 1 for baryons, and for 
a cosmological constant 7a = 0. 

The Einstein-Klein- Gordon equations that describe 
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this Universe are (in units c = h 

.2 / 4 



H = - 
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Pv 

Pb 



2 V 3 3 
3H* + V,* o = 0, 
4Hp z = 0, 
4H Pv = 0, 
3H Pb = Q, 



Pb 



with the Friedmann constraint 

.2 



H 



K 



(p*o + Pz + Pv + Pb + Pa) 



(4a) 

(4b) 
(4c) 
(4d) 
(4e) 



(5) 



being k 2 = 8irG, H = a/a the Hubble parameter and the 
commas stand for the derivative with respect to scalar 
field. Notice that background scalar quantities at zero 
order have the subscript 0. In order to solve the system 
of equations Q, we define the following dimensionless 
variables 



V6H' 
« \fpz 

K y/f>b 
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K vv 

K yfpl 

« Vpa 



Here we take the quadratic scalar potential V = 
™ 2 <i>g/2 with m ~ 1 x 10~ 22 eV. Using these variables, 
the equations Q for the evolution of the background 
Universe are transformed into 



—3 x — su + 




(7a) 
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sx + -n u, 




(7b) 






(7c) 






(7d) 


5<n-i) b, 




(7e) 
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Sq S , 




(7g) 



where the prime denotes a derivative with respect to the 
c-folding number N — In a, and n is defined as 



H 
IP 



(2x 2 



1 9 
~Z + 



-n. 



(8) 



We have introduced the variable s = Cq/H, where Cq 
is a constant, as a control parameter for the dynamics of 



< (see [57] for more details of the control parameter). 
With these variables, the density parameters f2i for each 
component i can be written as 

Q$ = x 2 + u 2 , 
^ z = z 2 , 



b 2 , 
I 2 . 



(9) 



subject to the Friedmann constraint 

x 2 + u 2 + z 2 + is 2 + b 2 + l 2 = 1. (10) 
In addition, we may write the EoS of the scalar field as 



fid 



(11) 



Since w$ is a function of time, if its temporal average 
tends to zero, this would imply that $ 2 -dark matter can 
be able to mimic the EoS for CDM. 

We solve the system of equations J7| for the back- 
ground Universe numerically with a four order Adams- 
Bashforth-Moulton (ABM) method. We take as initial 
conditions the best estimates from 5-years [51] and 7- 
( 6 ) years WMAP [52] values to ^ 0) = 0.73, = 0.22994, 



rt b 0) = 0.04, Oi 0) 



0.00004, n { ° } = 0.00002 and the ex- 
ponent k — 0. In Fig[l]we see the evolution of the scalar 
field $o- This figure shows how the scalar field oscil- 
lates very stark about the minimum of the scalar po- 
tential V, on the bottom of this figure we show a zoom 
of these oscilations. In Fig[2]we show the evolutions of 
the kinetic (<l> 2 /2) and potential energy (m 2 $ 2 /2) of the 
scalar field, as expected, the SF oscilations are translated 
into very stark oscilations for the kinetic and the poten- 
tial energies of the SF. However, observe the evolution of 
the dark matter density of the scalar field, that means 
/O* = ($ 2 +m 2 $ 2 )/2 in Figji] Note the following crucial 
point, although the kinetic and potential energies show 
very stark oscillations, the sum of both energies does not 
oscillates at all, this sum is the density parameter r2$ , 
which does not display any oscillation. The important 
point we have to take from now on into account is that 
the oscillations are not physical observables at all, they 
are a feature of the SF, what we observe in fact is the 
density of the SF which does not oscillates. 

Figj3] shows the numerical evolution of the density pa- 
rameters in our model 1 . At early times, radiation dom- 
inates the evolution of the Universe. Later on, the Uni- 
verse has an epoch where the energy density radiation is 



H. In Eq. (7g), sq is a constant and the exponent k is 



1 The scale factor is such that a 
the redshift zbya = (l + z)~ 1 . 



1 today, so that it relates to 



5 






K*„ 




FIG. 1: Evolution of the scalar field $o for the background 
Universe. The inset shows the fast oscillations of the scalar 
field. 



equal to the dark matter density, at z eq , then dark mat- 
ter begins to dominate the evolution. The recombination 
era in SFDM model occurs at z ~ 1000. At later times, 
the cosmological constant dominates the dynamics of the 
Universe at z\ ~ 0.5. The behavior is exactly the same as 
in the ACDM model. Furthermore, in the recombination 
era of the SFDM model, the neutrinos made up ~ 12 per 
cent of the Universe completely in agreement with the 
measurements of WMAP. Fig|4] shows the evolution of 
the EoS for the SF. Although the EoS varies with time 
(oscillations), the temporal average, (w$ ), drops to zero. 
Therefore, SFDM is like a pressureless fluid and behaves 
as CDM at cosmological scales ([55], [70] . [7Tj . [7J]). 



FIG. 2: Evolution of the kinetic (top panel) and potential 
(bottom panel) energy of the scalar field $o. 



Background Universe: The Fluid Approach 

In this section we perform a transformation in order 
to solve the Friedmann equations analytically with the 
approximation H << to. 

Here we take the scalar potential as V = to 2 <£> 2 /2h 2 + 
A<I> 4 /4, (we remain that we use units where c = H = 1 
unless it is stipulated). Then for the ultra-light boson 
particle we have that m ~ 10 -22 eV. 

Now we express the SF, $o ; m terms of the new vari- 
ables S and po, where S is constant in the background 
and po will be the energy density of the fluid also in the 
background. So, our background field is proposed as 

(ipo e~ imt l n + %e imt ' n ) 



$ 



where, 



(12) 



(13) 



and with this our SF in the background can be finally 
expressed as, 



$o = 2\/po cos^ - mt/h), 



(14) 



with this we obtain 



$o = A) 



Po 
Po 



cos(S — mt/K) 



2(S — m/ti) sin(jS — mt/K) 



(15) 



To simplify, observe that the uncertainty relation im- 
plies that mAt ~ K, and for the background in the non- 
relativistic case the relation S/m ~ is satisfied. Notice 
also that for the background we have that the density 
goes as (lnpo) = — 3-ff, but we also have that H ~ 10~ 33 
eV << to ~ 10 -22 eV, so with these considerations at 



hand for the background, in ( 15 ) we have 



$o = 4 ~^2~< 5 o sin 2 (5 - mt/K) 



(16) 



Finally, substituting this last equation and equation 
( 14 1 into ([!]) when taking A = 0, we obtain 



P<s> 



2 2 
77T, Jit 

2— polsin 2 ( S ~mt/h)+ cos 2 (S -mt/h)} = 2—p . 

h z h z 

(17) 




FIG. 3: Evolution of the density parameters f2, for the back- 
ground Universe. SFDM model mimics the standard ACDM 
behavior. 
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FIG. 5: Analytical evolution of the potential (top panel) and 
kinetic (bottom panel) energies of the scalar field dark matter. 
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shown in Fi g. |3] In what follows we will use the analytic 
expressions ( 18 ) and ( 19 1 instead of the numerical ones. 



FIG. 4: Evolution of the SFDM equation of state for the 
background Universe. 



Comparing this result with ^ we have that the iden- 
tity f2$ = 2m 2 p /ti 2 holds for the background, so com- 
paring with (17 1, 



SCALAR FIELD DARK MATTER 
FLUCTUATIONS 



/ 2/5o — sinfS 1 — mt/H) 
n 



/ Til 

\/2p — cos(S — mt/h). 



(18) 



(19) 



In terms of the two analytic results, we show the evo- 



lution of the kinetic and potential energies ( 18 1 and ( 19 ) 



in Fig. [5j where for the evolution we used the c-folding 
number N and the fact that a ~ t n — > t ~ e N ' n . 

Observe the excelent accordance with the numerical 
results in [71] for the kinetic and potential energies of 
the background respectively. Therefore, we expect that 
if we sum the analytical kinetic and the potential energies 
we will obtain the same behavior for the SFDM density 



Nowadays it is known that our Universe is not exactly 
isotropic and spatially homogeneous like the FLRW met- 
ric describes. There exist small deviations of this model 
that were the cosmological seeds to lead the large-scale 
structure formation in the Universe. In this section we 
compute the growth of the SFDM overdensities 5p<z> at 
the stage when the density contrast S = Sp®/ p§ a is much 
smaller than unity. It is believed that the Universe was 
almost uniform after inflation, with very small density 
contrast. As the Universe expanded, the small over- 
densities grew until they began to collapse, leading to 
the formation of structure in the Universe. Here we con- 
sider small enough deviations in the FLRW model, so 
that they can be treated by linear perturbation theory. 
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The Linear Regime of SFDM perturbations: 
Field approach 



The 



In the SFDM paradigm, the dark haloes of the galaxies 
are BECs composed of ultra light scalar particles with 
to ~ 10" 22 eV then it would be about 10 68 m GeV" 1 
scalar particles to follow in a single galaxy and the occu- 
pation numbers in galactic haloes are so big that the DM 
behaves as a coherent classical scalar field $ that obeys 
the Klein-Gordon equation (D 2 + m 2 //i 2 )$ = 0, where □ 
is the D'Alambertian. Therefore, to study the structure 
formation in the Universe we describe a non-interacting 
SFDM model such that: i) it can describe it more as a 
field than as particles and ii) we find a function that only 
depends on the three spatial coordinates and time. 

First, we introduce the perturbed metric tensor in the 
FLRW background, we only consider scalar perturba- 
tions. We then give the equation of energy-momentum 
conservation and the Einstein field equations for the per- 
turbed metric. 

By definition, a perturbation done in any quantity, is 
the difference between its value in some event in real 
space-time, and its corresponding value in the back- 
ground. So, for example for the SF we have 



where the background is only time dependent, while the 
perturbations also depend on the space coordinates. Sim- 
ilar cases apply for the metric; 



5oo = -a 2 (l + 20), 
goi = a 2 B,i , 

g tj = a 2 [(l-2ip)5 ij + 2E, ij ]. 



(21) 



Here ip is a perturbation associated to the curvature and 
E is asociated to the expansion. We will work under the 
Newtonian gauge, which is defined when B = E = 0. An 
advantage of using this gauge is that here the metric ten- 
sor g^ v is diagonal, and so the calculations become much 
easier. We will only work with scalar perturbations, vec- 
tor and tensor perturbations are eliminated from the be- 
ginning, so that only scalar perturbations are taken into 
account. Another advantage in using this gauge is that 4> 
will play the role of the gravitational potential, this will 
help us to have a simpler physical interpretation, i.e., 
both potentials <j> and i\) are then related. This metric 
has already been used in other works (see for example 
[9], [64] and [65]). 

The perturbed energy-momentum tensor for the SF 
can be written as the background value T = T it) plus 
a perturbation ST = ST(x tJ ') where x^ = [t,x l ], i. e. 



ST. 



(22) 



We now derive the perturbed evolution equations for 
the different quantities mentioned above; the scalar per- 



turbation and the scalar potential ip. For the per- 
turbed energy-momentum tensor, we get 

5T° = -Sp* = -($ <5$-$ 2 ^ + U,* Q 5$), (23a) 

ST? = -I (23b) 
a 

ST* = 5p$ = - - V, 9o S$)S). (23c) 



Where in equations ( 23 ) the dot denotes differentiating 



with respect to cosmological time t, which is related to 
conformal time by the simple relation 



— = a — , 
arj at 



(24) 



In the Newtonian gauge, the metric tensor g^ v becomes 
diagonal and from this, in the trace of the Einstein's 
equations, the scalar potentials ip and (f> are identical and, 
therefore, tp plays the role of the gravitational potential 



ip - 4> = o. 



(25) 



Usually this equation contains a term of anisotropic 
stress, which vanishes in the case of a SF. Altogether, 



(20) the perturbed Einstein's equations SG* = k 2 ST* to first 



order for a SF in the Newtonian gauge are 



— 8irGSpq> 
SirGSpq, 



!»//(. , + H(j>) - -^V 2 < 
or 



2(0 + 1/0),,, 

2[0 + 4ff0 + (2i7 + 3 J ff 2 )</>]. (26) 



These equations describe the evolution of the scalar per- 
turbations. For the evolution of the perturbations in the 
SF we use the perturbed Klein-Gordon equation 



5$ + 3HS$ - ^V 2 <5$ + V,^ 5$ + 2V,$ <t> - 4$ <; 



= 0. 

(27) 

In order to solve the equations (26) and (27 1, we will 



turn to Fourier's space. The beauty of this expansion re- 
laying on the fact that each Fourier mode will propagate 
independently. To first order, the derivation of Fourier's 
components is straightforward. The perturbation <5$ re- 
lates to its Fourier component S&k by 

6<P(t,x r ) = I cPkSQfak^expiikix'), 

d 3 M$ fc exp(ifc l x i ), (28) 



where k is the wave number. Here the wave number is 
defined as k = 2tt/L, being L the length scale of the 
perturbation. 



The perturbed equations ( 26 1 altogether with the SF 



transformation read 



Taking the time average of equation ( 35 ) we obtain 



87rG(3#<i>o<S$ fc ) 



2k 2 



2(H</> + 

2[4> + AH4>+ (2H + 3H 2 ) 



-87rG($ (5$ fe 

-<t>$l + V,« 5$ k ), 

(29a) 

87rG$ £$fc, (29b) 
87rG($ 5<i> fc 

_ ^*(J$ fc )(29c) 



5$+3H( < > - < w$ >)(5$ = 3</>< F$ > - < Grf, > . 
< dp$ > 

In equation ( 38 1 for the radiation and matter dominated 



eras the first term in the parenthesis goes as < 5p§ > / < 
<5p$ >ss 0, see for example [70], also we can see from Figj4] 
that (w$ ) — > 0. It is easy to see the temporal average of 



and the Klein-Gordon equation (27) transforms into 



5$ k + 3H5$ k + ^T<5*fc + V,^ (5$ fc - 4t6$ + 2</>V,$ = 0. 

(30) 

These set of equations describe the evolution of the per- 



the terms F<s> and G^ of the equation ( 38 1 . From Figj6] 
we observe that that (F$) tends to one. On the other 
hand, since we are using post-newtonian approximation 
(G^) drops to zero, meaning that the second term on the 
right-hand side of Eq. (38 1 dissapears. So, due to the 



scalar oscillations around the minimum of the scalar po- 
tential, the unperturbed SF behaves very similar to the 



turbations. Eq. (|29a|) makes reference to the evolution CDM model, but equation p) tell us that their pertur- 



of the energy density, Eq. ( 29b ) to the evolution of the 



gravitational potential and finally, Eq. ( 30 ) refers to the 
perturbations over the SF. 



Now, from the time derivative of (23a), we have 

<5p<t> = ($0 + Vi)5$fe 

+(5$k + V,^ 8$ k - $ o 0)$o ~ 20$ o $ o .(31) 

Performing a Fourier transformation of the above equa- 
tion and combining it with the unperturbed and per- 
turbed Klein-Gordon equations pb| a nd (30) respec- 
tively, and with the use of equation ( 29 ) we arrive at 



2k 2 

S Pi , = -6iJ<M$fc + H^lH - -^>{H4> + <j>) + 3^ 2 Q . 

(32) 

And on the other hand we have, 

Spt, + 5 P<S> = 2$ <5$ fc - 2§l4>, (33) 



Thus, 

5p<j> = — 3H (<5p$ + 8p^>) 



2k 2 



-{H(j> + <t>) +3<p®l (34) 



This last equation can be expressed in terms of the den- 
sity contrast 5$ making use of the equations from the 
background, we have 



, dp. 



5$ + 3H(^- — w$)<5$ = 30(1 + w$) — Gs, 
op® 



where we have defined the function G$ as 

2k 2 4> + H(f) 



a 2 n 2 p$ 



(35) 



(36) 



It is convenient to define the function F§ as well, 
F <j) = l+ uj$. 



bations too, all the growing behavior for the k modes are 
recovered and preserved so far. Finally, we conclude that 
Eq. ( 38 ) resembles the equation for the density contrast 
as in the standard CDM model [17] . This means that the 
SFDM perturbations in the linear regime grow exactly as 
CDM perturbations. 




FIG. 6: Evolution of (F<&) term involved on the right-hand 
side of equation ( 38 1 . 



The rescaled equations for SFDM perturbations 



We now study the cosmological evolution of our model 
through a dynamical system. In order to solve the system 
of equations for the perturbations given in Eqs. (29) 



(37) and (30), we define the following dimensionless variables 
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2;iven by 



U 



•i'2 



"-2 



A 



(39) 



Using the variables denned in Eq. ( 39 ) , the equations 



can be transformed into an autonomus system with re- 
spect to the e-folding number N, then the density con- 
trast can be written in terms of these variables 



6 = 



2[x(z2 — xl\) — Uzi] 



(40) 



Therefore the transition between the linear and non- 
linear regime of the SFDM perturbations can be studied 
through a numerical code. 



The Natural Cut of the mass power spectrum 



We rewrite equation ( 30 ) in the following form 



fc 2 

**fc+(-a +V,*« )<5$ fe = -3ff(5* fc +40*o-2<£V,* . (41) 



Equation ( 41 1 is the one of a harmonic oscillator with a 



damphing term 3H5<& and an extra force. This equation 
contains oscillating, growing and decaying solutions, de- 

7,2 

pending on the factor \ + If this factor is positive, 

we expect that the fluctuations oscillate, only if the fac- 
tor is negative we will have growing up fluctuations. This 
factor contains the scale fc 2 /a 2 , which is always positive. 
Thus, if we expect to have growing up fluctuation of the 
scalar field, this happen only when the potential V,$$ is 



negative. Then equation ( 41 1 will have growing up solu- 



tions for certain values of the scale k/a, this is, only for 
enough small values of the scale, which means for enough 
big fluctuations. This is the way we can fix the mass of 
the scalar field, the potential depends on the value of the 
scalar field mass. Observing the smallest structure we see 
in the Universe, we are able to fix the mass of the scalar 
field. Using this observations one can avoid the problem 
of substructure of the ACDM. In order to do this, in [70] 
the mass of the scalar field was fixed in 10 _22 eV. 



The Linear Regime of SFDM perturbations: The 
Fluid Approach 

The SF <E> has very stark oscillations from the begin- 
ning, these oscillations are transmitted to the fluctua- 
tions which apparently seems to grow very fast and are 



too big. Nevertheless, this behavior is not physical, be- 
cause we only see the oscillations of the fields, but we can 
not see clearly the evolution of its density, [72] . In order 
to drop out these oscillations, in what follows we perform 
two transformations. The first one changes the perturbed 
Klein-Gordon equation into a kind of 'Schrodinger' equa- 
tion and the second transforms this last equation into 
a hydrodynamical system, where we can interpret the 
physical quantities easier and the observable quantities 
become much clear. 

First, we consider the perturbed Klein-Gordon [27] 
where we have set <f> = 

<5$ + 3H5<S> - -r V 2 <5$ + 7,*$ 6® + 2V,* = 0. (42) 
a 2 - 

Now we express the perturbed SF <5$ in terms of the field 



,5$ = V]/ g-imt/H _|_ q,* e imt/h 



(43) 



term which oscillates with a frequency proportional to m 
and W = ^(x,t) which would be proportional to a wave 
function of an ensamble of particles in the condensate. 
With this equation and the expresion for the potential of 



the scalar field, (42 1 transforms into 



3 h 2 
HV) + — (□4- + 9A|4'| 2 *) + m</>1< 



2 ' 2m 
where we have defined 



□ 



dt 



lv 2 . 

r,2 



0, (44) 



(45) 



Notice that this last equation could represent a kind 
of ' Gross-Pit aevskii' equation in an expanding Universe. 



The only modification of equation ( 44 1 in comparison to 



the Schrodinger or the Gross-Pitaevskii equation is the 
scale factor a -1 associated to the co- moving spatial gra- 
dient and that the Laplacian V 2 = d\ transforms into 
the D' Alambertian □. 

To explore the hydrodynamical nature of bosonic dark 
matter, we will use a modified fluid approach. Then, to 
make the connection between the theory of the field and 
the condensates waves function, the field is proposed as 
(Madelung transformation [52"] 1 



* = 



<p e 



iS 



(46) 



where will be the condensates wave function with p = 
p(x,t) and S = S(x,t), [33]. Here we have separated 
\1/ into a real phase S and a real amplitude y/p and the 
condition 



VE"]/* = p is satisfied. From (46) we 



have 



hS jm 



iHp pas- 

m 



tup 

n 2 



2m 2 



VP 



-VSWp- 
2a 2 



m 

[V(Sym)] 



- —(s/my = 0. 



(47) 
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Now, taking the gradient of (47) then dividing by a and 
using the definition 



we have, 

h + 

V + 



-VS 



3Hp - — pUS + -vVp - — 'pS = 0, 
mam 



(48) 



Hv + J—;Vp+ -V<; 
2ap a 

1 



h 2 _ y ^dVp 



2m? a \ ^/ p 

-(n-V)ti - ftte + Hv){S/m) = 0. 
a 

where in (47) w = 9h 2 X/2m 2 and in (49) we have defined 



(49) 



p = wp'. 

It is worth noting that, to this moment this last set of 
equations do not involve any approximations with respect 



to equation ( 44 ) and can be used in linear and non- linear 
regimes. 

Now, neglecting squared terms, second order time 
derivatives and products of time derivatives in this last 
set of equations we get, 



dp 
dt 
dv 

at 



V-(pv) + 3Hp = 0, (50) 

Hv + (W)v - ^V(^V 2 p) + LoVp 

V0 = O, (51) 
4nGp, (52) 



where the equation for the gravitational field is given by 
the Poisson's equation (52). In these equations we have 
introduced r = a(t)x, such that 1/aV = V = <9 r . 



Equation (48) shows the proportionality between the 
gradient of the phase and the velocity of the fluid. Note 
that v can represent the velocity field for the fluid and 
p will be the particles density number within the fluid. 
Also there exists an extra term of third order for the 
partial derivatives in the waves amplitude which goes as 

this term would result in a sort 



the gradient of 2 ^ D ^ 



sfp 

of 'quantum pressure' that would act against gravity. We 
remain that <j) represents the gravitational field. These 



two sets of equations ( 50 ) and ( 51 1 would be analogous to 



the Euler's equations of classical 'fluids', with the main 
difference that there exists a 'quantum term', which we 
will call Q and will be given by Q = 2^^/=? which can 
describe a force or a sort of negative quantum pressure. 
For equation ( 50 1 we have that p will represent the 



mass density or the particles density number of the fluid, 
where all the particles would have the same mass. Fi- 
nally, these equations describe the dynamics of a great 
number of non-interacting identical particles that mani- 



fest themselves in the form of a fluid, also equation (44) 



can describe a great number of non-interacting but self- 
interacting identical particles in the way of a Bose gas, 



when the probability density is interpreted as the density 
number. 

Now, these hydrodynamical equations are a set of com- 
plicated non-linear differential equations. To solve them 
we will restrict ourselves to a vecinity of total equilib- 
rium. 

For this let po be the mass density of the fluid in equi- 
librium, the average velocity Vq will be taken as zero in 
equilibrium, so we will only have v(x, t) out of equilib- 
rium. Then, the matter in the Universe will be considered 
as a hydrodynamical fluid inside an Universe in expan- 
sion. This system will then evolve in this Universe and 
later on they will collapse because of their gravitational 
attraction. 



Then from ( 50 ) for the mass density of the fluid in 



equilibrium we have, 



dpo 
dt 



3H Po = 0, 



with solutions of the form 



Po 



Poi 



(53) 



(54) 



where as we know, in general if we have an equation of 
state of the form p — cop and consider CDM or dust 
as dark matter such that p = it holds that p cx a~ 3 . 
Then, when the scale factor was small, the densities were 
necessarily bigger. Now, the particles density number 
are inversely proportional to the volume, and must be 
proportional to a~ 3 , therefore the matter energy density 
will also be proportional to a -3 , result that is consistent 



with our expression (54). In addition, the cosmological 



behavior of po is in agreement with the numerical results 
show in Fig. [3] 

Now for the system out of equilibrium we have 



d5p 
~dt 
dSv 
~dt 
V 2 56 



3H5p + p W-Sv = 0, 



Hdv- 
4nG5p, 



1. 



V(^V 2 ^) + ujVSp + V54> = 0, 

(55) 



2m 2 ' '2 po 



equations that are valid in a Universe in expansion. In 
order to solve system ( 55 ) we look for solutions in the 



form of plane waves, for this the convenient ansatz goes 
as 

Sp = pi(t) exp( ik • x/a), 
dv = Vi(t) exp( ik • x/a), 
6<fi — 4>i{t) exp( ik • x/a). 

where x is the position vector and k is a real wavevector 
which corresponds to a wavelength L. If we substitute 
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these ansatz in the set of equations ( 55 ) , we then have 



dpi 
dt 
dvi 



+ 3Hpi 



i — k • V\ = 0, 
a 



at a \ pa 



4 ' G h 



0i + 47rGp-pi = 0, 
where we have defined the velocity 

, h 2 k 2 



(56) 
0,(57) 

(58) 
(59) 



To solve the system is convenient to rotate the coordinate 
system so that the propagation of the waves will be along 
the direction of one of the axes. For this we know that the 
velocity vector can be divided into longitudinal (parallel 
to fe) and transverse (perpendicular to fe) parts, such that 
V\ = vh + V2, where V2 is the vector perpendicular to 
the wave propagation vector fe • V2 = 0. In terms of V2 



for equations (56)-(58) we have 



dpi 
dt 
dv 
~dt 



3-ffpi 



Hv 



■ Po 



l — k V 

a 

,,2 



a po k z 



(60) 

■w)ft = 0, (61) 



in addition to an equation for V2, dv2/ dt+Hv2 = 0, with 
solutions V2 = C/a with C a constant of integration, i.e., 
perpendicular modes to the wave vector are eliminated 
with the expansion of the Universe. Now, if we use the 



result (54 1, then equation ( 60 ) can be written as 



_d_ 
dt 



ik 2 i> 



(62) 



System ( 60 )-( 61 ) can be treated as in the case of a Uni- 



verse with no expansion, so combining the two equations 



and with the aid of ( 62 ) , we get 
d 2 S 



dt 2 



2H- 



dS 
lit 



upo J =2 



AirGpo 



5 = 0, (63) 



where 5 — pi/po — Pi/Po is defined as the density con- 
trast. This will be a fundamental equation in the under- 
standing of the evolution of the primordial fluctuations. 



RESULTS 

First we will give a brief summary of the results for 
the ACDM model, this will enable us to make a direct 
comparison with our results. 

For CDM the equation for the evolution of the density 
contrast is given by, 



d 2 5 
~dl? 



2H- 



dS 
~dl 



,k 2 



AnGpo U = 



(64) 



FIG. 7: Evolution of the perturbations for the CDM model 
(dots) and SFDM model (lines) for k = 1 * 10" 3 /iMpc _1 . No- 
tice how after the epoch of equality (a eq ~ 10 -4 ) the evolution 
of both perturbations in nearly identical, a — 1 today. In this 
case we have taken A = 0. 



where c s is defined as the sound velocity (which in our 
case it is not). Now lets analyze equation (64) at the 



beginning of the matter dominated era a time just af- 
ter the epoch of equality, and just before recombination 
when the radiation has cooled down and the photons do 
not interact with the electrons anymore, for a relativis- 
tic treatment see |36| . In this era, a is a eq , practically 
all the interesting fluctuation modes are well within the 
horizon, and the evolution of the perturbations can be 
well described within the newtonian analysis. At this 
time, matter behaves like dust with zero pressure. So 
we have a ~ t 2 / 3 , c 2 k 2 /a 2 and p ~ t~ 2 therefore 
H = (2/3) l/t. For equation (64) we have 



6 = 0. 



cPS 4 1 dS 2 1 
d£2 + 3 1 "fa ~ 3^2- 

The solutions to this equation are of the form 

C2 
t ' 



5(t) -> t 2 ^C\ + 



(65) 



(66) 



where C\ and C% are integration constants, from this 
solution we can see that we have modes that will disap- 
pear as time goes by, and modes that grow proportionally 
to the expansion of the Universe. This is an important 
result, because then the density contrast will grow pro- 
portionally to the expansion of the Universe when this 
is dominated by matter. Then, these fluctuations can 
maybe grow and give life to the galaxies, clusters of galax- 
ies and all the large-scale structure we see now a days. 

Now lets see what happens to the SFDM at this epoch 
(a a eq ). The evolution of the perturbations in this case 



will be given by equation ( 63 ) 



In general we have that in equation ( 63 ) the term v q is 



very small throughout the evolution of the pertubations 



(v q < 10" 



for small fe), so it really does not have 



a significant contribution on its evolution. 
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FIG. 8: Evolution of the perturbations for the CDM model 
(dots) and SFDM model (lines) for k = 1 * 10~ 2 ftMf>c _1 and 
A / and negative. Notice how after the epoch of equality 



10" 



the evolution of both perturbations is now dif- 



ferent from the one in Fig. [jj a — 1 today. In this case we 
can clearly see that the SFDM fluctuations grow quicker than 
those for the CDM model. 



big in size, this means that they can only give birth to 
large structures. These fluctuations can then help for 
the formation of large clusters or other large-scale struc- 
ture in the Universe at its early stages (around a ^ a eq ). 
Then, as these kind of SFDM can only interact with ra- 
diation in a gravitational form it is not limited by its 
interaction with radiation, and the dark matter halos 
can then create potential wells that will collapse early 
in time giving enough time for the structures to form. 
Then if DM is some kind of SFDM, the luminous matter 
will follow the DM potentials giving birth to large-scale 
structure. 

Finally, a further analysis is required to obtain better 
results of the growth of scalar perturbations with self- 
interacting SFDM as well as to study the transition be- 
tween linear and no linear regime of scalar fluctuations. 
Nevertheless, we conclude that the SFDM model at cos- 
mological scale is a plausible alternative to DM nature 
because it preserves all the success of the ACDM model. 



When the condition A = is taken we can have a 
BEC that might be or might not be stable, if there exists 
stability the results of SFDM arc consistent with those 



obtained from CDM (in this case both equations ( 63 ) and 
(64) are almost equal), the condition of stability for the 
BEC in the SFDM case will come from the study of A 
together with Q. 



d 2 S 



2H- 



dS 

~dt 



^ a z 



6 = 0, 



(67) 



As we can see in Fig. [jj the perturbations used for 
the ACDM model grow in a similar way for the SFDM 
model, when A = 0, in this case both perturbations can 
give birth to structures quite similar in size, and this will 
happen with all the fluctuations as long as k is kept small. 

When A 7^ the results are quite different, so when dis- 
cussing the evolution of the density perturbations, there 
are two different cases: i) In the case of A > the ampli- 
tude of the density contrast tends to decrease as A grows 
bigger and bigger away from zero until the amplitude 
of the density takes negative values (around A ~ 10 8 ), 
telling us that this kind of fluctuations can not grow in 
time, and hence do not form a BEC. ii) On the other 
hand if A < the fluctuations for the density contrast al- 
way grow despite their size, this results means either than 
the fluctuations grow and form a stable BEC or than the 
density grows because it is collapsing into a single point 
and our BEC might be unstable, the study of the sta- 
bility of this fluctuations needs then to be studied with 
non-linear perturbation theory. These results are shown 
in Fig. [8] in both figures [7J and [8] the initial condition 
for S goes as S ~ 1 * 10~ 5 in accordance with the data 
obtained from WMAP. 

If these fluctuations result stable and because they are 



SCALAR FIELD DARK MATTER IN GALACTIC 
DYNAMICS 



In this section we study the SFDM/BEC model at 
galactic scales, we will focus on the cusp/core discrep- 
ancy between the numerical predictions of the standard 
model and the astronomical observations. 

As already said before, recent observations in far and 
nearby galaxies have shown that the ACDM model faces 
serious conflicts when trying to explain the galaxy forma- 
tion at small scales [32 ES]. ^ or instance, in the ACDM 
simulations the halos present rising densities towards the 



central region behaving as p 



well within 1 kpc 



[75] , On the other hand, several observations suggest 
that the dynamics (rotation curves) of dwarf and LSB 
galaxies are more consistent with a constant central den- 
sity [26l EH [89], this is most commonly known as the 
cusp/core problem 26 . 

Studying a wide range of galaxies of different mor- 
phologies and with magnitudes in the interval —22 < 
Mb < —8 Donato et al. [55] fit their rotation curves 
(RC) using a Burkert profile for the DM [TB] and found 
that 

\og(fj, a /M pc- 2 ) = 2.15 ± 0.2 (68) 
remains approximately constant, where 



Mo = Para 



(69) 



with po the central DM density and r$ the core radius. 
Similar result where found in Kormendy & Freeman [53] , 
Spano et al. [92]. Exploring further the constant value 
of fi for the DM, Gentile et al. [33] found that within 
ro the DM central surface density in terms of the mass 
inside it (M <ro ) is < £ >q,dm= M <r Jicr 2 w 72jJ£M Q 
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pc 2 ,the gravitational acceleration due to DM felt by a 
test particle at the radius ro was found to be 



gDM{r ) — Git < £ >o,dm= 3.2 



+ 1.8 
-1.2 



x 10 9 cm s 



(70) 

additionally they reported the acceleration due to the 



matter 

-2 



at r to be gw(?"o) = 5.7; 



-3.8 



luminous 
10- 10 cm s" 

In the ACDM model the galaxies have evolve through 
numerous mergers and grew in different environments, 
the star formation and basic properties of the galaxies 
are not expected to be a common factor among them. 
Therefore giving both the constancy of pg and the core 
in the central regions of galaxies seems very unlikely in 
this model. 

Any model trying to become a serious alternative to 
ACDM has to succeed in reproducing observations in 
which the standard model fails but also has to keep the 
solid description at large scale. It is always necesary to 
test the SFDM/BEC model with the two observations 
mentioned above, the cusp/core problem and the con- 
stant DM central surface density. In order to do this, 
in [85] it was used the Thomas-Fermi approximation and 
a static BEC DM halo to fit rotation curves of a set of 
galaxies. However, so far there was no comparison be- 
tween the density profile and the data, in that work this 
was done by fitting rotation curves of 13 high resolution 
low surface brightness (LSB) galaxies and additionally 
compare the fits to two characteristic density profiles 1) 
the cuspy Navarro-Frenk- White (NFW) profile that re- 
sults from N-body simulations using ACDM and 2) the 
Pseudo Isothermal (PI) core profile. The comparison al- 
low to show that the model is in general agreement with 
the data and with a core in the central region. For the 
second goal the meaning of a core is somewhat ambigu- 
ous. 

Large scale A^-body simulations using collissionless 
CDM the internal region of DM halos show a density dis- 
tribution described by a power law p ~ r a with a rj — 1, 
such behaviour is what is now called a cusp. On the 
other hand, observations mainly in dwarf and LSB galax- 
ies seem to prefer a central density going as p ~ r . This 
discrepancy between observation and the CDM model 
receives the name of cusp/core problem. Among the em- 
pirical profiles most frequently used to describe the con- 
stant density behavior in this galaxies are the PI [TU] , the 
isothermal [8j and the Burkert profile [16]. Even though 
their behavior is similar in the central region and is speci- 
fied by the central density fitting parameter, their second 
parameter called the core radius does not represent the 



same idea. For instance, in the PI profile (cq.( 79 )) we see 
that the core radius will be the distance in which the den- 
sity is half the central density. For the Burkert profile the 
core radius R b L urk will be when p bark (R burk ) = p b urk /4 
and for an isothermal profile (I) [§2] p I {R I c ) = /Oq/2 3 / 2 . 
Hence, we see an ambiguity in the meaning of the core 



radius, they get the same name but the interpretation 
depends on the profile. If we want to compare the cen- 
tral density of LSB galaxies with that of NFW, it usu- 
ally suffice to have a qualitative comparison, so far this 
is what we have been doing by fitting empirical profiles. 
However, high resolution rotation curves demand a more 
quantitative comparison. Indeed, if we want to test mod- 
els by fitting RCs we have to know the specific meaning 
and size of the core, then we will be able to tell if a model 
is consistent with a cusp or not by making a direct com- 
parison with the data. 

In |85j it was proposed a new definition for the core 
and core radius that allow us to decide when a den- 
sity profile is cusp or core. Using this definition in the 
SFDM/BEC model discussed above it is possible to find 
that the SFDM/BEC model can reproduced the constant 
value of po and as a crosscheck with the PI profile, in [55] 
was found that the results are in very good agreement 
with observations. This argues in favor of the model and 
the core definition. 



DM density profiles 

Following 85 we provide the dark matter density pro- 
files that will be used for the analysis. 

The case in which the dark matter is in the form of 
a static BEC and the number of DM particles in the 
ground state is very large was considered in Boehmer 
& Harko [14]. Following this paper and assuming the 
Thomas-Fermi approximation [HI [82] which neglects the 
anisotropic pressure terms that are relevant only in the 
boundary of the condensate, the system of equations de- 
scribing the static BEC in a gravitational potential V is 
given by 



= -pW, 



V 2 V = AnGp, 



with the following EoS 



p(p) = u oP 2 , 



(71) 



(72) 



(73) 



where Un 



2nh 2 a 



-, p is the mass density of the static 



BEC configuration and p is the pressure, as we are con- 
sidering zero temperature p is not a thermal pressure but 
instead it is produced by the strong repulsive interaction 
between the ground state bosons. Assuming spherical 
symmetry and denoting R as the radius at which the 
pressure and density are zero, the density profile takes 
the form [Ml 



Pb{t) 



B sin(fcr) 



(74) 
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where k = ^/GmFJWa = tt/R and p^ = ps(0) is the 
BEC central density, m is the mass of the DM particle 
and a is the scattering lenght. The mass at the radius r 
is given by 



m(r) 



47rpo f sin(fer) 



A- 2 



kr 



cos(fcr 



(75) 



from here the tangential velocity V of a test particle at 
a distance r, is 



T/ 2 (r) = 



47rGpjf / sin(fcr) 



cos(fcr 



fc 2 y fcr 

The logaritmic slope of a density profile is defined as 

d(log p) 



(76) 



d(log r) 



using (74 1 in ( 77 1 we obtain 



a{r) 





' 7rr\ 


l--cot( 







(77) 



(78) 



We study now the Pseudo Isothermal (PI) profile. The 
empirical core profiles that exist in the literature fit two 
parameters, a scale radius and a scale density. A charac- 
teristic profile of this type is 



Ppi 



n PI 
Po 



1 



(r/Rc 



(79) 



this is the PI profile [TU] . Here R c is the scale radius and 
Pq 1 is the central density. The rotation curve is 



V(r)pi = jAnGp^Rl ( 1 - — arctan 



(80) 



The NFW profile emerges from numerical simulations 
that use only CDM and are based on the ACDM model 
[3"0l 1751 177] . In addition to this, we have chosen this 
profile because it is representative of what is called the 
cuspy behavior (a « —1) in the center of galaxies due to 
DM. The NFW density profile [77] and the rotation curve 
are given respectively by 



PNF\v(r) 



P< 



(r/Rs)(l + r/R s y 



(81) 



V N Fw(r) = y/AnGpiRh 



In 1 



r 

R. 



r/R s 



1 



r/Rs 
(82) 

Pi is related with the density of the universe at the mo- 
ment the halo collapsed and R 2 S is a characteristic radius. 

To solve the core radius ambiguity and to unify the 
concept for future comparison, in [55] it was found that 
a good definition for the core is a region where the density 
profile presents logarithmic slopes a > — 1 and the core 



radius will be the radius at which the core begins, that 
is to say, for radius smaller than the core radius we will 
have a > — 1, this means that its value r' is determined 
by the equation 



a{R) = -1. 



(83) 



The advantages of this definition are that the interpreta- 
tion is independent of the profile chosen (also notice that 
it applies to the total density profile and is not restricted 
to that of DM) and in virtue of the same definition we 
can directly tell if a DM model profile is cored or cuspy. 
With this new definition the specific distance at which 
the core radius occurs still depends on the profile chosen 
but now the physical interpretation is only one. In the 
following when we refer to both the core and core radius 
we adopt the previous interpretation. 



Applying the definition to ( 74 ) we get the core radius 



for the SFDM/BEC profile Rb and for comparison we 



use ( 79 ) because in turns out that the parameter R c cor- 
responds to the core radius as defined above. Finally, 
fitting the NFW profile provides a direct comparison be- 
tween a cusp and core and hence to the cusp/core prob- 
lem. 



SFDM vs. PI and NFW 



We see from (|74| that the SFDM/BEC model satis- 
fies p ~ r° near the origin, but a priori this does not 
imply consistency with observed RCs. Therefore we re- 
produce the fit of the profiles to thirteen high resolution 
observed RCs of a sample of LSB galaxies given in [55] . 
The RCs were taken from a subsample of de Blok et al. 
[27] . we chose galaxies that have at least 3 values within 
~ 1 kpc, not presenting bulbs and the quality in the RC 
in Ha is good as defined in McGaugh S. S. |73j . The RCs 
in this work omit galaxies presenting high asymmetries 
and included in the error bars are experimental errors in 
the velocity measurement, inclination and small asymme- 
tries. Because the DM is the dominant mass component 
for these galaxies we adopt the minimum disk hypothesis 
(neglects baryon contribution to the observed RC). 

As the difference between a core and a cusp is not over- 
lapped only for data values inside lkpc and given that in 
the interval ~ 1 to lOkpc the slopes of core and cusp pro- 
files are very similar, which can lead to the wrong con- 
clusion that cuspy halos are consistent with observations, 
we determined the logarithmic slope and the uncertainty 
following de Blok et al. [27] with the difference that we 
fit only the data within 1 kpc and that there is no need 
of an uncertain "break radius" . 

In Fig. [9] we show the fits to the RC data and the den- 
sity profiles, also shown are the core radius in the BEC 
(magenta) and in PI (blue) profiles. The gray arrow is 
the fit that determines a, the length denotes the fitted 
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region and is bounded by R\ that denotes the nearest 
radius to lkpc where a data point is given. We show 
a by fitting values inside R\ and we also show the core 
radius for the SFDM/BEC profile Rb in order to com- 
pare it with R c (see [55] for more details of the fitting 
parameters). 



Discussion 

The fits of the RCs in Fig [9] prove that the solid-body 
like behavior characterized by a linear increase in the ve- 
locity in the central region is more consistent with the 
core PI and SFDM/BEC profiles than the cuspy NFW. 
The fits within Ri give an average value of 0.27 ± 
0.18 consistent with those obtained in de Blok et al. [27 
a=-0.2 ± 0.2 and with a=-0.29 ± 0.07 reported by 
Oh Se-Heon et al. [72] analyzing 7 THINGS dwarf galax- 
ies. The case of ESO1870510 might be considered to be 
consistent with NFW profile, however it is the innermost 
value that considerably decreases a, being an irregular 
galaxy more central data near the innermost region is re- 
quired to discard the possibility of any violent event that 
might have caused such a slope value. 

The density profiles corresponding to the RCs fits are 
also shown in Fig [9] for each galaxy. We see that the 
SFDM/BEC fit slightly deviates from the farthest data 
points as a result of having a finite radius R that is fixed 
by the same data. This discrepancy is due to the fact 
that the halo might be more extended than the value R. 
As a matter of a fact, the more the extended the "flat" 
outer region in the RCs the more conspicuous the dis- 
crepancy. Harko & Madarassy [35J recently found finite 



temperature corrections to ( 79 1 , this suffice to alleviate 



the latter problem in LSB galaxies and dwarfs but not 
for bigger galaxies. Some solutions have been proposed 
including vortex lattices [HH I104j and adding more nodes 
g7J ET] in the solution of system ([FT] and ([72])), but so 
far no final conclusion has been reached. 

When comparing the SFDM/BEC and PI core radius 
we find a general difference of ~ 2kpc, the core size in the 
PI is approximately 2 kpc smaller than the BEC core size, 
but the PI central density is larger. In U4115, U11557 
and U11583 both profiles are very similar which results 
in a similar core and central density values, this can also 
be taken as a consistency check for our core definition. 

Comparing the values of i?s we did not find a tendency 
to a common value. Assuming that the core radius deter- 
mines the transition where the DM distribution changes 
from the outer region to the inner constant central den- 
sity, the lack of the unique value means that there is not 
a common radius at which this transition takes place. 



For the second test [83 uses Rb to calculate (69). We 



have already seen that R c and Rb are generally differ- 
ent and Rb is not a fit parameter. Hence a priori Rb 
is not expected to correlate with With the fitting 



parameters it was obtained 

log(/zf /M QP c~ 2 ) = logpf R B = 2.05 ± 0.56, (84) 



log(/^ J /M qP c-') = log R c = 2.08 ± 0.46, (85) 
for the average values in the SFDM/BEC and PI pro- 



files respectively. We see the excellent agreement of ( 84 ) 



with ( 85 1 that was used as a crosscheck and with ( 68 1 in 



which a much bigger sample was used. The agreement 
has shown that the SFDM/BEC model is capable of re- 
producing the constancy of the value po, something that 
because of the cuspy nature is not possible in the NFW 
profile. 

In Fig. [10] we plot the above values for each galaxy. If 
we define the DM central surface density mentioned in 
the introduction, we have for the BEC profile 



< s >0,DM~ M <R. 



(86) 



where M < r b is obtained from (75) evaluated at Rb- We 
show that for U11748 the value logp B is considerably 
above the rest and with the largest uncertainty. For this 
reason, in the following analysis we omit both, this value 
and the smallest one that corresponds to ESO1200211. 
Doing this we get < S > B DM ~ 191.35 M /pc 2 , and 
for the acceleration felt by a test particle located in Rb 
due to DM only we have <7dm(-Rb) ~ 5.2 xl0~ 9 cms -2 
broadly consistent with ( |70"| . 

The fact that all galaxies present approximately the 
same order of magnitude in qdm{Rb) might suggest that 
Rb represents more than a transition towards a constant 
density, it can give us information about the close rela- 
tion between DM and the baryons. Moreover, in view 
of the lack of a unique core radius, we can interpret the 
transition in DM distribution as an effect of crossing a 
certain acceleration scale instead of a radial length scale. 
Such interpretation reminds us that given in MOND but 
with the big difference that the acceleration scale found 
is for DM and is not a postulate of the model. 

To determine which interpretation causes the transi- 
tion whether an acceleration scale or a length scale we 
will need to study the properties of larger samples of 
galaxies with the new telescopes. 



CONCLUSIONS 

In this work, we revisit an alternative paradigm to dark 
matter nature known as scalar field dark matter or Bosc- 
Einstcin condensate dark matter model. In this model a 
fundamental scalar field plays the role of dark matter in 
the Universe. The hypothesis is that this scalar field un- 
dergoes a phase transition very early in the Universe lead- 
ing to the formation of Bose-Einstein condensate drops. 
Therefore, the dark haloes of the galaxies in the Uni- 
verse are huge drops of scalar field. To explore if this 
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FIG. 9: Observed LSB galaxy rotation curves and density profiles with the best halo fits. Below each RC is its density 
profile along with the fits. Shown are the PI (green dashed-line); the BEC (red solid line) and NFW (black dobble-dotted 
line) DM halo profiles, the observational data is drawn with error bars. The gray arrow denotes the best fit to the data 
within Ri and the vertical arrows denote the PI (blue) and BEC (magenta) core radius. 
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FIG. 10: Plot of log (/ij? /M e pc~ 2 ) and log (/^ /Mqpc' 2 ) 
for each galaxy. N is an identifier for each galaxy. Here we 
observe that these values remain approximately constant in 
both profiles, this also serves as a crosscheck for our defini- 
tion of R B in the SFDM/BEC profile. The green dashed-line 
represents the mean values given in (841 and ( 85 1. 
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hypothesis is a viable alternative to the standard model 
we developed both cosmological and galactic studies. 

In the cosmological regime, we have shown, with a field 
approach, that the SFDM/BEC model with an ultralight 
mass of 10 -22 eV mimics the behavior of the cosmological 
expansion rate predicted with the ACDM model. More- 
over, although in general a scalar field is not a fluid, it 
can be treated as if it behaved like one. Thus, we per- 
formed a fluid approach for the cosmological evolution of 
the scalar field and we found that the analytical expres- 
sions for the kinetic and potential energies of the scalar 
field is in excellent agreement with our previous numeri- 
cal results. 

The interesting cosmological behavior of the scalar 
field indicates that their scalar fluctuations can be the ap- 
propriate for the purpose of structure formation, because 
overdense regions of SFDM/BEC can support the forma- 
tion of galactic structure. Thus, we have revisited the 
growth of SFDM/BEC fluctuations in the linear regime. 
Within the linear theory of scalar perturbations we ob- 
tain an equation for the evolution of the density contrast. 
This equation differs from the density contrast equation 
for CDM, however the extra terms tend to the values 
of the standard equation. Therefore, as a goal of this 
model, we obtained that the scalar perturbations grow 
up exactly as the CDM paradigm. In addition, we study 



the liner growth of the scalar fluctuations with a fluid 
approach. In the case of the SFDM/BEC with A — we 
have shown that for the matter dominated era and for big 
structure this model simulates the behavior of CDM be- 
cause in general in a matter dominated Universe for low- 
k, v q tends to be a very small quantity tending to zero, 



so from ( 63 ) we can see that on this era we will have the 



CDM profile given by (|64j), i.e., the SFDM density con- 
trast profile is very similar to that of the ACDM model 
(see Fig. [7]). On the contrary for A ^ both models have 
different behavior as we can see from Fig. [8] The numer- 
ical results suggest that linear fluctuations on the SFDM 
can grow faster than those for CDM around a ~ 1CU 2 . 
Here an important point is that although CDM can grow 
it does so in a hierarchical way, while from Fig. [8] we can 
see that SFDM can have bigger fluctuations just before 
the ACDM model does, i.e., it might be that no hierar- 
chical model of structure formation is needed for SFDM 
and is expected that for the non-linear fluctuations the 
behaviour will be quite the same as soon as the scalar 
field condensates, in a very early epoch when the energy 
of the Universe was about ~ TeV. Thus, the standard 
and the SFDM/BEC model can be contrasted in their 
predictions concerning the formation of the first galax- 
ies. If in the future we see more and more well formed 
and massive galaxies at high redshifts, this could be also 
a new indication in favour of the SFDM/BEC paradigm. 

On the other hand, we also studied the implications of 
a SFDM/BEC model at galactic scales. We find that the 
SFDM/BEC model gives a constant density profile that 
is consistent with RCs of dark matter dominated galaxies. 
The profile is as good as one of the most frequently used 
empirical core profiles but with the advantage of com- 
ing from a solid theoretical frame. We fit data within 
lkpc and found a logarithmic slope a = — 0.27± 0.18 in 
perfect agreement with a core. It is important to notice 
that the cusp in the central regions is not a prediction 
that comes from first principles in the CDM model, it 
is a property that is derived by fitting simulations that 
use only DM. In addition, we review a new definition for 
the core and core radius that allows a definite distinction 
when a density profile is core or cusp. Using this defini- 
tion we find the core radius in the SFDM/BEC profile to 
be in most cases over 2kpc bigger than the core radius in 
the PI profile. As a second result of the core definition, 
we show the constant value of /irj which is proportional to 
the central surface density. This result is one of several 
conflicts that jeopardize the current standard cosmolog- 
ical model. 

Finally, with all these intriguing results SFDM/BEC 
model could be a serious alternative to the dark matter 
problem in the Universe. Although the observational ev- 
idence seems to be in favor of some kind of cold dark 
matter, if we continue to observe even more galaxies at 
higher redshifts and if higher resolution observations of 
nearby galaxies exhibit a core density profile, this model 
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can be a good alternative to ACDM. We expect that the 
close future observations in galaxies surveys can decide 
the nature of the dark matter. 
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